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ABSTRACT 

Active galactic nuclei (AGNs) are well-known to exhibit flux variability across a wide range of 
wavelength regimes, but the precise origin of the variability at different wavelengths remains unclear. 
To investigate the relatively unexplored near-IR variability of the most luminous AGNs, we conduct a 
search for variability using well sampled JHK S -band light curves from the 2MASS survey calibration 
fields. Our sample includes 27 known quasars with an average of 924 epochs of observation over 
three years, as well as one spectroscopically confirmed blazar (SDSSJ14584479+3720215) with 1972 
epochs of data. This is the best-sampled NIR photometric blazar light curve to date, and it exhibits 
correlated, stochastic variability that we characterize with continuous auto-regressive moving average 
(CARMA) models. None of the other 26 known quasars had detectable variability in the 2MASS 
bands above the photometric uncertainty. A blind search of the 2MASS calibration field light curves 
for AGN candidates based on fitting CARMA(1,0) models (damped-random walk) uncovered only 7 
candidates. All 7 were young stellar objects within the p Ophiuchus star forming region, five with 
previous X-ray detections. A significant y-ray detection (5a) for the known blazar using 4.5 years of 
Fermi photon data is also found. We suggest that strong NIR variability of blazars, such as seen for 
SDSSJ14584479+3720215, can be used as an efficient method of identifying previously-unidentified 
y-ray blazars, with low contamination from other AGN. 

Subject headings: quasars: general surveys: 2MASS 


1. INTRODUCTION 

The temporal flux variability from active galactic nu¬ 
clei (AGN), detectable in nearly all wavelength regimes, 
contains information on the underlying emission pro¬ 
cesses a nd sou r ce geome try that is otherwise difficult to 
probe (IUlrich et al.lll99?t) . However, precise details of the 
physical mechanism generating the observed n ucle ar vari¬ 
ability in AGN remain unclear (lAntonuccil 12013T ). Cur¬ 
rent and future large-scale photometric time-domain sur¬ 
veys have motivated many recent studies of the optical 
broadband variability properties of various AGN sub¬ 
classes using large numbers of well-sampled light curves. 
This has been especially useful for AGN identification 
and selection dKellv et al .1120091 iKozlowski et al.l^OlObt 
IZu et al.ll20T3h. 

Beyond the optical, large-scale surveys of AGN 
variability have been pursu ed at many other wave¬ 
lengths, incl uding the radio (iThvaearaian et all l201lh . 
ultraviolet (iGezari et all I2013D . and y-ray regimes 
(lAckermann et aLl[2l) 1 2 i. Fewer studies of AGN variabil¬ 
ity have focused on the infrared (IR), due in part to the 
expectation that non-variable dust and the host galaxy 
dominate the emissio n at these wavelengths (e.g. see 
iKishimoto et al.ll2008f) . Explanations of AGN variability 
in the optical as due to localized tem pera t ure fluctua tions 
in inhomogeneous accretion di sks dR uan et a U 120141 1 or 
global accretion rate changes dPerevra et al.l l2006f) pre¬ 
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diet little flux variability from disk emission in the IR. 
Reprocessing of variab le optical emission by a dusty torus 
has been suggested bv lSuganuma et a l. (2006) as a possi¬ 
ble cause for the NIR variability, although the time-lags 
between the optical and NIR a re difficult to detect. For 
example iMcHardv et al.l (|2007f ) have shown small time 
lags between NIR and X-ray data for the AGN 3C273, 
with the NIR variability leading the X-ray by ~1.5 days, 
implying this NIR variability originates near the nucleus. 

A study of NIR variability i n the Spitzer Dee p Wide- 
field Survey Bootes field by IKozlowski et, al.l (l2010alf 
showed that although only 1.1% of objects appear to be 
variable in the near-IR (NIR), the vast majority of the 
variable objects were AGN. Small samples of AGN moni¬ 
tored in both the optical and NIR have shown that while 
nuclear variability is prevalent, the amplitude of varia¬ 
tions decreases towards l onger wavelengths. For exam¬ 
ple, iHonig fe Kishimotol (120111) have produced realistic 
models of the optical to NIR variability for NGC 4151, 
which show a small amplitude time lag and decreasing 
amplitude variability as a function of wavelength. How¬ 
ever, this is not true for blazars, whose NIR flux is dom¬ 
inated by non-thermal emission from a relativistic jet, 
rather than circum-nuc lear dust whic h is not expected 
to be non-variable dCutri et al.1 119851) . Recent results 
from optical/NIR monitoring of Fermi 7-ray blazars 
have shown that flat-spectrum radio qu asars are more 
variable in the NIR than in th e optical (Bonning et_al.| 
l2012t iSandrinelli et al.l I2013D . ISandrinelli et al.l (120141 ) 
have used 7 year photometric light curves of 7 blazars, 
spanning the optical to NIR bands ( VRIJHK ), to find 
that variability amplitudes increase with wavelength for 
these objects. 

In this study we search for NIR variability of luminous 
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Fig. 1.— WISE color-color space for the Cal-PSWDB sources 
with WISE detections (black points and grey filled contours), 
known AGN (red cro sses), and t he W ISE blazar strip (black 
open contours) from IChoi et~ . I|2014D . The known blazar 
SDSSJ14584479+3720215 is highlighted (blue circle). 

AGN using well sampled light curves from the Two Mi¬ 
cron All Sky Survey (2MASS) calibration data, finding 
one highly variable blazar (SDSSJ14584479+3720215). 
We explore the utility of NIR variability-based selection 
of AGN, including its use in identifying counterparts to 
Fermi 7 -ray sources. We also study the nature of non- 
AGN contaminants with similar NIR colors and variabil¬ 
ity properties, laying groundwork for the AGN variability 
science possible with current and future multi-epoch IR 
surveys. 

2. DATA AND AGN SELECTION 

The 2MASS survey observed the full sky in the near-IR 
over the timespan of Jun e 199 7 to February 20 01 using 
the J, H , and K s bands dSkrutskie et al.ll2006D . Photo¬ 
metric zero-points for calibrating the survey were based 
on hourly observations of a set of 40 standard Helds, each 
815 x 1° in size, and spaced evenly throughout the sky. 
These fields were visited repeatedly over the course of the 
survey, resulting in 562 to 3692 epochs of observation per 
field. This produced JHK S light curves for 113,030 indi¬ 
vidual objects, known as the 2MASS Calibration Point 
Source Working Database (h ereafter Cal-PSWDB. iCutril 
120061: IPTavchan et al.ll2008bl ). These light curves are the 
most precise NIR time domain survey to date and pro¬ 
vide a novel dataset in which to study the infrared vari¬ 
ability properties for a wealth of astropliysical phenom¬ 
ena. This dataset has produ ced the best samp led NIR 
light curve for a RR Lyr star (|Szabo et al.ll2014D. a hunt 
for s tellar flares at long wavelengths dDavenport et al.l 
20121) . a large sample of young stellar objects (YSO’s, 
Plavchan et, al.ll2008al iParks et al.ll2014D. and many well 

characterized bina ry star systems dBecker et al.l 120081 
I Quillen et al.l 120141 Davenport 2014 in prep). Here we 
make use of these data to search for NIR variability of 
AGN with unprecedented precision. 

The Wide-field Infrared Survey Explorer (WISE, 
IWright et al.ll2010h has created the deepest mid-IR sur¬ 
vey of the entire sky in four photometric bands, from 
3.4 gin to 22 gm. Several recent studies have demon¬ 
strated the utility of WISE colors for effectively separat¬ 
ing AGN from stars and unresolved galaxies. We spa¬ 


Fig. 2.— J-band (top) and J — K s color (bottom) light curves 
for the known blazar, SDSSJ14584479+3720215. The median pho¬ 
tometric error is shown for comparison (vertical bar). 

tially matched the Cal-PSWDB objects to the WISE all 
sky data release using a match radius of 1"5. This re¬ 
sulted in 52,148 objects with JHK S light curves and a 
detection in at least one WISE filter. The distribution of 
these point sources in WISE color-color space is shown 
in Figure [Q For co mparison , we s how the WISE “blazar 
strip”, defined in I Choi et al.l (120141) as the Gaussian Ker¬ 
nel Dens it y Es ti mate of the WISE colors of blazars from 
iMassaro et al.l (120091) . 

To find previously known AGN in our data, we spa¬ 
tially matched these 52K objects to the “Milliquas Cat¬ 
alog,” version 3.E0. We recovered 27 known AGN with 
light curves in Cal-PSWDB, 25 of which had WISE col¬ 
ors shown in Figure [I] These AGN had between 26 
and 3482 epochs of Cal-PSWDB photometry, with an 
average of 924 epochs. Of these, 19 have been spec- 
troscopi cally confirmed b y the Sloan Digita l Sky Survey 
(SPS S; lYork et all 120001 : IBovv et all 120111: IPa.ris et all 
f 20 ll) . one of which has been identified as a blazar 
(SDSSJ14584479+3720215). 

3. AGN VARIABILITY IN THE NIR 

This set of known AGN with well sampled Cal-PSWDB 
light curves creates a unique dataset in which to search 
for AGN variability in the NIR. To determine if these 
objects were variable in the 2MASS bands, we computed 
the error-weighted root mean square variability for ea ch 
light curve, E, following the prescription in iSesar et al.1 
(|2007l) . Objects were classified as variable if they had 
Ej > 10 -2 mag. Of the “Milliquas” AGN in our dataset, 
only the known blazar, SDSSJ14584479+3720215, exhib¬ 
ited variability above this threshold. The J-band and 
J—K s color light curves for this blazar, with 1972 epochs 
of Cal-PSWDB photometry spanning over three years, 
are presented in Figure [5] The remaining 26 known AGN 
in our sample displayed no significant NIR variability us¬ 
ing this metric. 

We found strong variability in 

SDSSJ14584479+3720215 in all three 2MASS bands 
with our densely sampled light curves. However, as 
seen in Figure [2] no significant variation in J — K s color 

5 http://quasars.org/milliquas.htm 
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was found for this object, indicating the underlying 
variability was “grey” over this wavelength regime. A 
larger sample of such well studied variable objects is 
needed to test if this is truly a generic property of 
blazars in the NIR. 

Characterizing photometric variability timescales has 
proven to be a very efficient means of selecting 
a clean sample of AGN from time domain sur¬ 
veys (e.g. IKozlowski et alJl2010bt iMacLeod et alJ l2011t 
iButler fc Blooml 120111 ). The so-called “damped random 
walk” (DRW) model provides a robust fit to AGN light 
curves as a function of two observable free parameters: 
the characteristic timescale r, and the variability ampli¬ 
tude SFqo. Stars are easily removed when fit with this 
model, as their primary sources of variability are typi¬ 
cally stochastic (e.g. flares or flickering) or strongly peri¬ 
odic (such as pulsations or rotation) in nature. Com¬ 
bining DRW variability models with multi-wavelength 
photometry provides even gr eater p ower in separating 
blazars from “normal” AGN dRuan et alil2012l ). 

From our Cal-PSWDB-WISE matched sample we 
selected all light curves having at least 50 good J- 
band epochs (PH_QUAL=A,B or C), median J-band 
magnitudes brighter than 16.5 mag, and WISE colors 
of W 1 — W 2 > 0.3 to remove main sequence stars 
(iDavenoort et al.ll2014H . We also required sources to be 
detected as variable using the Xj > 0.01 cut described 
above, and removed objects known to be periodic (bi¬ 
naries and_gulsating ^ variables) or long period variable 
stars (P lavchan et, alJ [2008al Davenport 2014 in prep). 
This yielded 22 variable Cal-PSWDB targets to search 
for DRW timescales, which included our known blazar 
object, but did not select any other previously identified 
AGN. For completeness we also analyzed the other 26 
known “Milliquas” sources for DRW variability. 

For each object we down-sampled the light curves to 
a single data point per night, using the median of the 
J-band photometry, which reduced the light curves to 
an average of 144 epochs. We then ran t he DRW fitting 
code f rom IKozlowski et alJ (l2010bfl and IMacLeod et alJ 
(2010) on the reduced J-band light curves for all 47 ob¬ 
jects. We used the p robability thresholds outlined in 
IMacLeod et alJ (1201(1) to determine which light curves 
had significant DRW variability. The known blazar had 
a weakly constrained DRW timescale of logr = 2.0j^'g 
days, and a significant driving amplitude of log a = 
— O.lljto’^ 9 nrag year -1 / 2 . The relative likelihood for 
the DRW versus an infinite or unconstrained timescale 
was \og(Luke/^inf) =1.4, where values greater than 
0.05 are considered significant. Similarly, the relative 
likelihood for the DRW model versus pure noise was 
log(L;ifc e /Lnoise) = 153, which was much greater than 
the threshold for significance of 0. None of the other 26 
previously identified AGN showed any signs of significant 
DRW variability. However, 7 of the variability-selected 
objects did show significant DRW variability. We discuss 
the origin of these objects in the following section. 

We additionally fit these light curves using the Contin¬ 
uous Auto-Reg ressive M oving Average, or CARMA(p,q), 
models of iKellv et ail ( 2014 ). where p is the auto¬ 
regressive order and q the moving average order. The 
DRW is a special case of the CARMA models, equiv¬ 
alent to a CARMA(1,0), or a first-order continuous 


auto-regressive process with no moving average. In the 
CARMA formalism, a DRW model’s power spectral den¬ 
sity is described as a single Lorentzian function with 
a characteristic break frequency. As mentioned above, 
the DRW has been shown to empirically model the be¬ 
havior of QSO light curves well, but is otherwise in¬ 
flexible. By looking at higher order CARMA models 
we search for more complex behavior in the light curve 
power spectral density, such as stochastic variability on 
multiple timescales or quasi-periodic oscillations. For 
SDSSJ14584479+3720215 we found a moderate prefer¬ 
ence for the CARMA(2,0) model compared to a DRW 
when using the deviance information criterion for model 
selection. The DRW model provides a sufficient fit to the 
data, in that the sequence of residuals is independently 
and normally distributed. A Markov Chain Monte Carlo 
analysis of the CARMA(1,0) model provided 1-er confi¬ 
dence levels on the timescale of31<r<115 days, and 
on the standard deviation of the driving noise of between 
0.14 and 0.20 magnitudes. 

4. YSO CONTAMINATION 

The recovered DRW timescales for the seven 
variability-selected objects we found in the previous sec¬ 
tion had an average timescale of logr = 1.05 days, sig¬ 
nificantly shorte r than the typical tim escale seen in the 
optical for AGN (IMacLeod et alJl2010D , as well as for our 
blazar target. This was reproduced in the CARMA(1,0) 
models as well. These seven objects all reside in field # 
90009, which was centered on the star forming region p 
Ophiuchus, making these likely YSOs. Of these YSOs, 
five had X-ray detections in the literature, one coming 
from a targeted X- ray monitoring campaign of the field 
(jGagne et alJ 120041) . and four fro m the Chandra X-ray 
source catalog fE vans et~aPl2010H . The other two YSOs 
in our sample did not have a published X-ray counter¬ 
part within 0.1°. NIR light curve properties for YSOs 
within this field in the Cal-PSWDB hav e be en studied in 
great detail by iPlavchan et all (l2008al) and iParks et~akl 


These seven YSOs had mid-IR colors that placed them 
within the blazar strip, NIR variability amplitudes com¬ 
parable to our known blazar, and light curves best pa¬ 
rameterized by a CARMA(1,0) model. We thus highlight 
YSOs as a potentially important source of contamination 
for classifying AGN from future multi-wavelength time 
domain searches. While using longer wavelength filters 
can impro ve the efficiency of se lecting AGN from normal 
stars (e.g. iMassaro et al l 120121) . YSOs remain a poten¬ 
tial contaminant f or object class ifica tion (see also color 
spaces in Fig 7 of iKoenig et al.ll2012f) . However, these 
objects may be distinguished from bona fide AGN based 
on their ~10X shorter characteristic timescales under a 
CARMA(1,0) model. YSOs are also mostly confined to 
the Galactic disk and dense star forming regions, and as 
such can be avoided by surveys. 

5. A 7 -RAY DETECTION OF SDSSJ14584479+3720215 

The Fermi Gamma-ray Space Telescope provides the 
deepest survey to date in the 100 MeV to 100 GeV 
regime, and has disc overed 886 y-ray AGN in the 2nd 
Fermi AGN Catalog (lAckermann et afll2011l) . The vast 
majority of these AGN are blazars. However, ~30% of 
Fermi sources in the Fermi 2nd Point Source Catalog 
















































4 


Davenport et al. 


200Mev - lOOGeV Counts Map 


Fitted Model 



232 230 228 226 224 222 220 218 
RA [deg] 


232 230 228 226 224 222 220 218 
RA [deg] 


II 


i 


60 

54 

48 cT 
42 01 

30 § 
24 El 
18 b 

l/l 

12 

6 

0 


QJ 36 
O 


Difference Image 




232 230 228 226 224 222 220 
RA [deg] 


CD 


1.2 

0.8 

0.4 "D 
in 

0.0 § 
o 
u 

0.4 

0.8 

- 1.2 


Fig. 3.— Left: 200 MeV - 100 Gev counts maps of ~4.5 years of Fermi data, centered on the known blazar (circled). Middle: Simulated 
counts map for the best fit model of all sources in the field, including the known blazar, detected with TS = 30. Right: Difference image 
between the observed and model counts maps. 


12FGL: iNolan et ahll2012f) are unidentified, largely due 
to the poor angular resolution of Fermi s pair-conv ersion 
Large Area Telescope (LAT; I Atwood et al.1 [20091. All 
reliably-associated AGN in the 2nd Fermi AGN Cata¬ 
log are confidently detected by Fermi with Test Statistic 
TS > 25, corresponding to approximately > 5tr detec¬ 
tions. For a Fermi counts map of a Region of Interest 
(ROI) centered on a new possible source, the TS detec¬ 
tion significance is the defined to be the ratio of like¬ 
lihoods TS = 2(logL(source) - logL(nosource)), where 
L (source) is the likelihood of a model of all known sources 
in the ROI and a new source at the central location, 
while L(nosource) is the likelihood for the same model 
but no new source at the central location (|Nolan et al.1 
120121) . The detection significance is approximately y/TS. 
Thousands of potential faint 7 -ray blazars lie below the 
TS > 25 (> 5<t) detection threshold, but are difficult 
to localize due to extremely low photon counts. These 
faint 7 -ray blazars have large positional uncertainties, 
with error ellipses with radius >20' for 3cr Ferm i detec¬ 
tions using 4 years of photon data (B allet et al.1 120131 : 
iThompson et al.ll2014D . 

The strong JHK S variability exhibited by the blazar 
SDSSJ14584479+3720215 relative to the quasars in our 
sample suggests that this NIR variability is likely to be 
indicative of strongly beamed jets, which are character¬ 
istic of blazars, and are known to produce strong 7 -ray 
emission. This allows for NIR-variability based identi¬ 
fication of previously unidentified or undetected 7 -ray 
blazars. SDSSJ14584479+3720215 was not included in 
the Fermi 2FGL catalog, and thus had TS < 25 is 2 
years of photon data. However, motivated by the strong 
NIR variability we observe, we investigated whether this 
blazar can be detected at high significance using more 
photon data. We downloaded approximately 4.5 years 
of available Pass 7 Fermi photon data of the Source 
event class, spanning Mission Elapsed Time 239557417 
to 383570220 seconds. We selected photons in the 200 
MeV to 100 GeV energy range, within a 10° Region of 
Interest (ROI) centered on the blazar. 

Using the standard binned likelihood analysis scripts 
provided in the Fermi Science Tools package, we per¬ 
formed photon event selection with a ROI-based zenith 
angle cut of <100° and a rock angle cut of <52°. We 
produced a spatial model of the likely 7 -ray sources 
within 15° of the blazar based on the 2FGL catalog, 
using the 2FGL best fit values for the template spec¬ 


tral model of these sources. We included the blazar 
SDSSJ14584479+3720215 as an additional point source, 
modeled as a simple power-law spectrum. The spec¬ 
tral normalization for all sources within 15° was al¬ 
lowed to vary, and all spectral parameters were set as 
free parameters within the ROI. We included the ap¬ 
propriate 2-year Pass 7 Galactic diffuse emission and 
extragalactic isotropic diffuse emission models in our 
source model fit to the observed counts map. The blazar 
SDSSJ14584479+3720215 is detected with TS = 30, 
approximately a 5<r detection, with a 7 -ray power-law 
spectral slope of 2.40 ± 0.18, typical of Fermi blazars 
(|Ackermann et ahll2011f) . This detection significance is 
consistent with the TS ~ 30 found in a preliminary 
version of the Fe rmi 4-year point source catalog by 
I Ballet et al. (120131 b produced using 4 years of photon 
data. We note minor discrepancies in our modeling of 
the brightest object in the field, noticeable in the lower- 
left corner of the difference image. These are well-known 
issues in Fermi source analysis caused by the finite resolu¬ 
tion of the counts maps, and do not affect our detection. 

6 . DISCUSSION 

We have conducted a search for NIR variability from 
AGN using the 2MASS C al-PSWDB light curves. No 
significant variability was found for 26 of the 27 known 
AGN in our sample. One previously known AGN, showed 
strong variations in the 2MASS filters, and we have pre¬ 
sented the 1972 epoch NIR photometric light curve for 
this blazar (SDSSJ14584479+3720215). This previously 
known blazar is a benchmark object, displaying almost 1 
magnitude of variability in the JF[K S bands, but with no 
significant variation in J—K s color. This is the best sam¬ 
pled NIR light curve of a blazar ever measured. We note 
a handful of other AGN candidates from the 2MASS Ca l- 
PSWDB were recently identified bv 1 Quillen et al.1 d2014li . 
including the blazar SDSSJ14584479+3720215. 

Our detection of SDSSJ14584479+3720215 as a 7 -ray 
source, selected by infrared variability, has implications 
on the association and classification of multi-wavelength 
counterparts of Fermi sour ces. C urre nt selecti on meth¬ 
ods such as rnid-IR color dMassaro et a 1.1 .2012 : and op¬ 
tical variability (|Ruan et al.l l2012h suffer from low effi¬ 
ciency. This is often due to “normal” Type 1 quasars, 
which have similar colors in the mid-IR and are also vari¬ 
able in the optical. As this pilot study demonstrates, 
these issues may be surmounted by the inclusion of IR 
variability information, since normal AGN are not sig- 
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nificantly variable at this wavelength regime. 

A search for other such light curves in the 2MASS Cal- 
PSWDB reveals a large number of contaminating YSOs. 
This study establishes a baseline for the AGN variability 
science possible with current and future multi-epoch IR 
missions such as the WIS E, the Wide Field I nfrared Sur¬ 
vey Telescope (WFIRST; iGreen et al.l I2012T), and Vista 
Variables in the Via Lactea ('VVV: iMinniti et aLll2010i ) 
surveys. 
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TABLE 1 

NIR Properties of the 27 “Milliquas” AGN found in the Cal-PSWDB. Our benchmark blazar is indicated (f). 


2MASS ObjectID 

Nepoch(^) 

(J) 

{J - K s ) 

W1 

1+1 - 1+2 

1+2-1+3 

J002410.86 - 015646.6 

1063 

16.86 

1.45 

14.91 

1.33 

3.55 

J015454.88 + 004044.6 

30 

16.96 

1.43 

15.80 

1.41 

3.05 

J015451.41 + 005933.2 

339 

17.00 

1.51 

15.29 

0.42 

2.71 

J015429.74 + 002711.0 

86 

17.03 

1.52 

15.61 

1.27 

3.38 

J034105.84 + 070917.4 

11 

17.08 

1.18 

15.39 

1.48 

2.39 

J055708.14 - 002414.3 

254 

16.99 

1.48 




J085116.85 + 120028.6 

2580 

16.72 

1.75 

13.58 

0.97 

2.85 

J094235.96 +591354.7 

68 

17.08 

1.49 

15.22 

1.49 

3.27 

J094249.43 + 593206.6 

778 

16.23 

1.56 

13.39 

0.87 

2.58 

J122125.95 - 001006.1 

1873 

16.68 

1.41 

14.89 

0.51 

2.91 

J122131.10 + 000742.7 

9 

17.01 


15.20 

1.32 

3.20 

J122144.68 - 001141.5 

19 

17.09 

1.76 

15.78 

1.43 

3.44 

,1122152.13 + 001719.0 

49 

17.02 

1.41 

15.77 

0.59 

3.18 

,1121427.03 + 350907.9 

2737 

16.29 

1.10 

13.38 

1.24 

2.85 

,1121418.95 + 352920.6 

18 

17.10 


16.50 

1.20 

3.84 

,1121408.54 + 355021.5 

1212 

16.91 

1.56 

14.63 

0.20 

2.64 

.1144057.24 - 000951.1 

916 

16.40 

1.13 

13.80 

1.59 

2.97 

,1150026.22 - 005428.2 

1568 

16.79 

1.64 

14.16 

1.18 

3.09 

.1150043.39 - 005820.7 

1500 

16.50 

1.48 

13.88 

0.74 

2.92 

.1145846.09 + 371235.9 

58 

17.06 

1.38 

15.32 

0.98 

3.56 

,1145844.84 +372021.8t 

1927 

16.33 

1.77 

14.03 

1.10 

2.85 

,1163124.44 + 295301.8 

1191 

15.51 

1.73 

12.46 

1.07 

2.52 

J204110.28 - 052626.3 

21 

16.99 

1.49 




J220037.69 + 211051.9 

888 

16.75 

1.32 

14.34 

1.45 

2.90 

J220028.41 +203902.1 

15 

16.88 


15.22 

1.20 

2.36 

J231808.32 + 001152.3 

7 

16.99 

1.46 

15.06 

1.39 

2.80 

J231821.13 + 002937.2 

358 

16.95 

1.52 

14.62 

0.86 

2.41 






